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Abstract 

We study directed polymers subject to a quenched random potential in d 
transversal dimensions. This system is closely related to the Kardar-Parisi- 
Zhang equation of nonlinear stochastic growth. By a careful analysis of the 
perturbation theory we show that physical quantities develop singular behav- 
ior for d — > 4. For example, the universal finite size amplitude of the free 
energy at the roughening transition is proportional to \/4 — This shows 
that the dimension d = 4 plays a special role for this system and points 
towards d = 4 as the upper critical dimension of the Kardar-Parisi-Zhang 
problem. 
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I. INTRODUCTION 



The field of non-equilibrium growth processes has attracted quite a lot of interest in the 
recent past The simplest nonlinear model describes a growing surface on sufficiently large 
length scales as a height profile h{r, t) over a d-dimensional reference plane parameterized 
by r. The dynamics of this surface is given by the Kardar-Parisi-Zhang (KPZ) equation 

^^h = uAh + ^\{Vhy + r] (1) 

with a Gaussian white noise defined by 

T]{r, t) = and r/(r, t)r/(r', t') = a^5'^{r - r')5{t - t'). (2) 

The surface described by (|^) can be in different phases which depend on the dimension-less 
coupling constant 

and the space dimension d. In less than two dimensions there are two different phases: 
the weak coupling phase for = and the strong coupling phase for g ^ 0. Above two 
dimensions, there exists a critical value gc- The surface is in the weak coupling phase for 
\g\ < gc and in the strong coupling phase for \g\ > g^. Precisely at \g\ = gc the system 
undergoes a roughening transition. Whereas the linear growth equation for g = can easily 
be solved in any dimension, it is much more difficult to get information on the critical 
behavior of the other phases. 

The morphology of the surface in the different phases is characterized by the asymptotic 
scaling of the height-correlation-function 

mn,h) - Kr2,h)f) ~ In -r2p^/(t|ri -rsr^), 

which defines the roughening exponent x ^^id the dynamic exponent z. For g ^ ^ these 
exponents satisfy the hyperscaling relation z + x = 2 

In the strong coupling phase it is known that z = 3/2in(i=l |^ and z = 2 in the limit 
of infinite dimension In general dimensions all exact methods fail and only numerical 
and mode-coupling results are available, but they become less reliable in higher dimensions. 



Therefore there is still a very controversial discussion [§-|T0| about the existence of a finite 
upper critical dimension of the KPZ problem, i.e. a dimension above which the dynamical 
exponent z has the constant value 2. 

In the language of the renormalization group, the different phases belong to different 
fixed points. For less than two dimensions, there is one unstable fixed point at (? = 0, which 
governs the weak coupling phase, and one stable fixed point at — — oo, which governs the 
strong coupling phase. Above two dimensions, the weak coupling fixed point also becomes 
stable and a new fixed point describing the roughening transition appears in between the 
other fixed points. 



It has been shown that the strong-coupling fixed point is inaccessible by a perturba- 



tion expansion around (7 = 0. The situation is somewhat better for the fixed point describing 
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the roughening transition. The singularities, which arise in the perturbation series above 
two dimensions, can be treated in a systematic expansion [12,13| with parameter 



(4) 

In the framework of this e-expansion, one finds the exponents z* = 2 and x* = 0, which are 
exact to all orders in perturbation theory |]lT|. Hence, hyperscaling is preserved. 



However, this perturbation expansion breaks down for (i — > 4 since new singularities in 
the coefficients of the series expansions arise at e' = 0, where 

e'^6+l = i^. (5) 

The treatment of these singularities is the aim of this paper. We develop a systematic 
way to extract the behavior of physical quantities as — > 4 from the divergent series in 
e'. We show that in contrast to the singularities in e, these singularities translate into a 
non-analytic behavior of observable quantities. This stresses the importance of = 4 for the 
KPZ problem and therefore shows that four dimensions is a good candidate for an upper 
critical dimension. 

We will address this question not in the KPZ picture but by using the exact Hopf-Cole- 
mapping [0,11^ of the KPZ problem to a directed polymer in a random medium. A directed 
polymer (also called a string) is a line with a preferred direction, which is governed by its 
line tension. The energy of a given conformation is 

n[r] = \r\t) + Az/-3/2r/(r(t), t)dt, (6) 

where r{t) is a (i- dimensional vector which denotes the elongation of the directed polymer 
at a position along the preferred ("time") axis t, L\\ is the projected length of the directed 
polymer and rj is the random potential which appears in the KPZ equation. 

To study the structure of the new singularities in the perturbation series, we will use an 
especially simple physical quantity. We nevertheless expect, that more complicated quan- 
tities such as correlation functions show the same type of singularities. We will study a 
system where the projected length L\\ of the directed polymer is infinite while its transver- 
sal "motion" is restricted to a finite volume of a characteristic width L±_ = L}^"^ . Within 
this system we will concentrate on the dimension-less averaged free energy per unit "time" , 
which we call 



C((?,L)= hm —F{g,L,L\\). (7) 

The infrared regularization by the length scale L has to be introduced, because the series 
expansion starts at the Gaussian fixed point {g = Q), which has no intrinsic length scale. L 
moreover serves as the fiow parameter of the renormalization group considered below. 

Since at the roughening transition hyperscaling is preserved at least around 2 dimen- 
sions , therefore there are no corrections to the 1 / L\ = 1/L behavior of the free energy 



per length F{g,L,L\\)/L\\, and the finite size amplitude 
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C{g)= \iu,C{g,L) (8) 



L 



is a universal quantity. Near d = 2, one finds [|TT| for its value C* at the unbinding transition 

C*(e) ~e + 0(e2). (9) 

C is not only one of the simplest physical quantities to be calculated in our system, but 
also one of the most fundamental ones: it plays a role very similar to the central charge in 
two-dimensional conformal models |T7 |. 



The directed polymer problem with randomness described by can be treated analyt- 
ically via the replica trick. This means that it can be expressed as the limiting case N —>■ 
of directed polymers with a short-ranged interaction potential (see e.g. |T6|): 

•^0 ^ i=l i<j 

By d^), the potential is always attractive {g < 0) and captures the short-ranged correlations 
of the disorder (H). 

The finite size amplitude (|^) of the free energy per unit "time" of the random system 
can then be obtained from the limit of the finite size amplitudes per unit "time" 

and per pair of directed polymers 

Cig,L,N) 



A^(A^- 1) 



X hm ^[F{g,L,N,L\\) -F(0, L, iV, L,,)] 

of the A^-polymer system without randomness. 

The development of a new regularization scheme for the perturbation theory of C near 4 
dimensions is done in two steps. In the first step, the fact is used that for only two directed 
polymers independent techniques exist, that allow an exact solution of the 2-polymer- 
problem. Those methods (the transfer matrix approach and the resummation of the per- 
turbation series) are reviewed and extended to our system with a transversally restricted 
"movement" in section |I| of this paper 

In order to be able to get results about perturbation series where no exact solution exists 
any more, a regularization scheme is developed by comparison with the exact solution. The 
regularization scheme itself then is independent of the existence of the exact solution in the 
sense that it extracts the results directly from the coefficients of the new singularities in 
the perturbation series. It is shown that only the most divergent terms in every order of 
perturbation theory contribute to the first order terms of the limiting /3 function. 

As the second step this regularization scheme is then applied to the case of an arbi- 
trary number of directed polymers in section [TTTI , which can only be solved exactly in one 
dimension [jl8| . A diagrammatic expansion of the partition function and of the free energy is 



developed. The series expansion is calculated up to the fifth order in the coupling constant 
and many simplifications appear, that suggest a simple underlying structure. 

Since the most important diagrams with respect to the regularization scheme turn out 
to be just 2-polymer- diagrams, the main properties of the /? function can still be extracted 
from the general structure of the perturbation series. 
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Independent of the number of directed polymers in less then 4 dimensions there are two 
zeroes of the (3 function: one is the finite size amplitude C* of the free energy per unit "time" 
at the unbinding transition; the other one is at the Gaussian fixed point. In 4 dimensions the 
j3 function becomes just a straight line with a fixed slope. C* and the zero at the Gaussian 
fixed point coincide there. We find that near 4 dimensions the physical quantity C* behaves 
non-analytically as 

C*{e') ~ Vi'. (10) 

Since the behavior is independent of the number of directed polymers, it should also 
describe the scenario in the replica limit of a vanishing number of directed polymers and 
therefore the behavior of the weak-coupling fixed point of the KPZ problem. 

In section the main results are summarized again. Many of the technical details have 
been postponed to various appendices to keep the route of argumentation straight. 



II. TWO DIRECTED POLYMERS 

Since the problem of two directed polymers (A^ = 2) is especially simple, we start exam- 
ining this problem and try to generalize our results in the next section. 

The reason for the simplicity of the two polymer problem is that one can separate the 
problem in the motion of a "center of mass" (ri(t) + r2(t))/2, which just fluctuates freely 
independent of the coupling constant g and the problem in the relative coordinate f = 
r2 — Fi. Since the KPZ-problem is equivalent to directed polymers with a relative short- 
range interaction (\^(ri,r2) = g5'^{ri — r2)), we have to deal in the relative coordinate with 
one directed polymer, which is interacting with a straight line defect at the origin of the 
perpendicular space. The transition probability from r' to r within the projected contour 
length t (the restricted partition function) is given by the single path integral 



^t(r',r) 



/•(t,r) 




/ Vr exp 


[-L 


/(0,r') 





'''-'r\t') + V{r{t'))dt' 



(11) 



The problem of one directed polymer interacting with a defect line can now be treated 
in different manners: On the one hand the path integral formulation of the directed poly- 
mer problem is equivalent to a "quantum mechanics" and can be solved in analogy to the 
corresponding quantum mechanical problem. This is called the transfer matrix approach 
because the Hamiltonian of the quantum mechanics is the continuum transfer matrix of the 
problem. 

The second approach just sets up a perturbative expansion of the problem with the defect 
with the strength of the defect as the expansion parameter, which can either be calculated 
coefficient by coefficient or in the special case of two directed polymers can even be summed 
up completely. 

A comparison of these approaches shows that they lead to the same results. Because we 
know the exact solutions, we can use them to develop a treatment of the perturbation series 
in the interesting case e' ^ 0, which we afterwards can apply to problems where no exact 
solution is available any more, as for example the case of more than two directed polymers. 
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A. Transfer matrix results 

Since the path integral (|TT]) is formally a quantum mechanical path integral in imaginary 
time, the partition function obeys the Schroedinger equation 

(-A + K(r))Z,(r',r) = -|z,(r',r), (12) 

of a particle in a potential V{r), which can be attacked by standard quantum mechanical 
methods. 

If the movement is restricted to a finite volume characterized by some length L± = , 
the long-time (large projected length of the directed polymer) behavior of the partition 
function is an exponential decay, the decay rate of which is given by the ground state energy 
of the quantum mechanical problem 

(-A + \/(r))^(r) = EoV^(r), (13) 

The free energy per unit "time" for long polymers is therefore given by the ground state 
energy plus exponentially small corrections. The finite size coefficient can be extracted by 
studying the behavior of the free energy for large system sizes L_|_. At the Gaussian fixed 
point we expect from dimensional analysis that i?o ~ = The prefactor of this decay is 

C. In order to calculate C, we extend here the calculations in [|l^] to the case of "movement" 
in a finite volume. 



1. Precise definition of the model 

In order to solve the model in arbitrary dimensions, we have to use potentials and 
boundary conditions, which have a spherical symmetry. In this case, we can separate the 
wave function in a radial and an angular part. As it is well known from quantum mechanics, 
it is convenient to define the radial function $(r) to be r^'^^^^^'^ times the radial part of the 
total wave function. This function obeys for the lowest (zero) angular momentum - as it 
should be the case for the ground state - the radial Schroedinger equation 

- 1^*^^^) ^ (^-3)^(c^-l) J_^^^^ +V{r)^{r)=E^^{r). (14) 

We implement the attractive potential at the origin as a well potential with a small but 
finite extension a: 

VW = (-yf<" (15) 

a < r 

By defining 0(r) = <l>(ra), Eq = o^Eq, V{r) = a^V^ra) and Vq = o^Vq, the ultraviolet cutoff 
a can be eliminated and everything is written in dimension-less variables. 

The boundary conditions imposed to the directed polymer have to be radially symmetric, 
too. The two possibilities studied here in detail are on the one hand the Dirichlet bound- 
ary conditions, where the directed polymer is put into a spherical box of radius Lj_ with 



6 



hard walls, which means (f){L±/a) = 0, and on the other hand the von- Neumann boundary 
conditions, which impose that the first derivative of the radial part of the wave function van- 
ishes at the boundary of a spherical box of radius Lj_, in order to mimic a kind of periodic 
boundary conditions. This is equivalent to the condition 



L d_ 

a dy 



m = ) (16) 

y a 



for the radial function (p. 

Moreover we will consider harmonic boundary conditions. That means that we allow a 
"movement" in infinite space, but add an harmonic well potential (7/2)r^, which effectively 
restricts the motion. The associated width of the box is from dimensional analysis y/j = 
{7t/2){1/L\). (The factor of 7r/2 has only been introduced for convenience). 



2. Results for the free energy 

For the potential as given by |T^ and the first two boundary conditions given in the above 
section, it is relatively easy, to calculate the ground state energy, because the potential is 
piecewise constant. Only standard techniques of quantum mechanics are applied. A sketch 
of these calculations is given in appendix ^ 

Here, we just give the results for the ground state energy. Let us first consider the 
Dirichlet boundary conditions. 

For the free system (Vq = 0) the ground state energy is exactly given by 



a 



Eo = xf, [- ) , (17) 

with x\e\ being the smallest positive root of the Bessel function J\^\. This is also the asymp- 
totic behavior for < Vq < K- At the phase transition point (Vq = K) we have 



yf.i if:)' 



lei < 1 



4(|e|-l)(^) " |e|>l 

where y\e\ is the smallest positive root of J-\e\- 

For von-Neumann boundary conditions the ground state energy is zero at Vq = 0. At 
< Vo < K the asymptotic behavior is 

4|e|(|e| + l) ^a^2(H+l)_ ^^^^ 



^1 I J|.|(v^) _i\L 



VoJ|.|+i(VW) 

The decay is faster than quadratic which means that the coefficient we are looking for remains 
zero. This scaling behavior of the free energy per length in a system of finite transversal 
size at the unbinding transition is consistent with the scaling behavior of the eigenenergy of 
the bound state in the infinite system . 
At the phase transition we get 
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En 



VH(N^-i)(it) 



|e| < 1 
|e| > 1 



(20) 



where z\^\ is the smallest positive root of The decay for |e| > 1 is again faster than 

quadratic. 

The differences between the decay coefficients at the phase transition point and in the 
free case (which we defined to be C* in the introduction) are plotted against the dimension 
in Fig. 0. 
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FIG. 1. Transfer matrix results for the renormalized coupling constant. The dependence on 
the dimensionality is shown for von-Neumann (a) and Dirichlet (b) boundary conditions. 



It can be seen that C* stays finite as d approaches 4 with Dirichlet boundary conditions, 
whereas for von-Neumann boundary conditions C* approaches zero. From the first three 
terms of the series expansion of the Bessel function (the third term is the first one 

with a non- vanishing coefficient in the limit |e| — > 1) we get near |e| = 1 a behavior of 
z\,\ ~ 4y^(l + |e|)(l - |e|) +4(|e| - 1), which yields equation (0). 

For harmonic boundary conditions the differential equation in the outer area is a bit 
more complicated. Its solution is 

^y) ~ y--2W ,(V7aV) (21) 

with the Whittaker function W\^^{z). The condition of differentiability at ?/ = 1 is also more 
complicated to evaluate than for the von-Neumann or Dirichlet boundary conditions, but in 
the end we get 

at < < K and 
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at the phase transition Vq = V^. Therefore 

C = n{2-d) = 2ne. 



(24) 



B. General perturbation theory 

Starting from the partition 

H = + ^ '^it)dt with = 5'^{r{t)) (25) 

Jo 



into the free part (kinetic part and infrared regularization) and an interaction part we can 
easily set up a perturbation series of the partition function as 

A = (e-^/o'"*Wd*)^ (26) 

^0 



l + E(-^7)"/ dti.-.y dt„($(ti)...<l>(t„,))c 



^ 0<ii<...<t„<L| 

where the expectation values (. . .)o are taken with respect to Hq and Zq is the partition 
function just for Tig- Here the factor of n! has been canceled out by introducing a time 
ordering. 

The free energy is the negative logarithm of the partition function and can also be 
expressed as a series in g by systematically expanding the logarithm. Since all correlation 
functions are translationally invariant in the limit of large L\\, one of the integrations just 
gives a factor of L\\ which cancels out since we are interested in the free energy per length. 
As a result of this series expansion we get to the first orders 

L||-+00 /;|| 

,2 



g{m)o - g' / {mm)o - (O)) dt + 0{g'). 
Jo 

Since the structure of the singularities in the case d ^ 4 does not show up in the second 
order, we have to extend the perturbation series to higher orders. (Higher terms of this 
formula can be found in appendix p.) 

For two directed polymers with a S interaction it is now particularly simple to calculate 
the time ordered correlation functions, because they can be reduced to two-point-functions. 
Consider the propagator Gt{ri, r2) of the relative coordinate of the free (only infrared regu- 
larized) directed polymers. With this, the multi-point-function is expressed as 

, Idn!dr,GtAru0)Gt,-tA0,0)---Gt^-t.-A0,0)GL,-tMr2) 

/dri/dr2GL||(ri,r2) 

For long chain length we can drop end effects and assume ti and L\\ — 1„ to be large enough 
so that we can combine those terms to a one point function 
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. . . $(t„))o = Gt.^t, (0, 0) . . . (0, 0) ($(0))o (28) 

If we here set n = 2 we get the relation between the return probabihty Gt{0, 0) and the two 
point function, so that we can express the multi-point function by two point functions as 

X(j.(ii)t(i2))„...{t((„_i)4(i„))j. 

Using this formula, we can express all the coefficients of the perturbation series as mul- 
tiple integrals over products of two-point functions. It is even more convenient to introduce 
the dimension-less connected two point function as 

_ moMLs))o,c _ ($(0)$(Ls))o-($(0))2 
9i[s) = 



If we now introduce a dimension-less coupling constant 

uo = gL' = gLmo))o, (29) 
the finite size coefficient of the free energy per unit "time" reads 

C(„„, = L lim £Mf£M . (30) 



Uo-Uq gi{s)ds 
Jo 

+ul / giisi)gi{s2) - giisi + S2) dsids2 
Jo Jo 



+0(4), 

where the time-ordered integration has been broken up into integrals over the respective 
time differences. 

As we can see, the finite size coefficient is equal to the dimension-less coupling constant 
Uq to first order. Singularities in this series can only arise in the integrals of the higher order 
terms. The renormalized coupling constant must be defined such that all singularities in 
the higher order terms are canceled. This is done by imposing the renormalization point 
condition for the renormalized coupling constant 

ur{uo) = C{uo). (31) 

This renormalization point condition has the advantage, that it gives a physical meaning to 
the renormalized coupling constant. 

Now we have to study the perturbation series (^) in detail. Already in the third order 
the essential difference between the regularization scheme at e = and at e' = becomes 
obvious. Near e = 0, the second term in the third order coefficient does not diverge at 
all, so that only the first term, which obviously factorizes to the square of the second order 
coefficient, produces divergences. This statement is true to all orders of perturbation theory, 
which produces just a geometric series of divergences. The especially simple structure of 
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those divergences guarantees that the (3 function calculated to the second order is exact to 



all orders of perturbation theory ||TT|j21 



Near e' = the situation is totally different, because here also the second term diverges. 
So in every order combinations of different types of divergences occur, which do not factorize 
any more. 

Since the terms in the perturbation series are becoming quite nasty in higher orders, it 
is convenient to introduce a graphical representation. Therefore we will draw in the n-th 
order n points in a line which represent the n time stamps involved. For each function gi we 
then draw a line connecting the two time stamps which are the arguments of the two point 
function. With this representation we can sketch our series expansion up to higher orders 
as it is done in appendix 0. 

Careful inspection shows that all diagrams appear that have at least one line passing 
each interval between two points and do not have two lines leaving one point in the same 
direction. The corresponding prefactors up to the seventh order are given by 

(-1)*^-'-^ n v-^ (32) 

peP 

where / is the number of lines in the diagram, P is the set of points in the diagram and for 
each p E P Ip denotes the number of lines which pass the given point. 

This perturbation series can be further simplified by applying some relations among the 
integrals over two point functions as for example 



Those relations can be produced by a general mechanism via Laplace transformation as 
explained in detail in appendix y. 

Applying them consequently enables us to get rid of all "nested" diagrams. Since we are 
finally interested in the inverse series to this one, we can just invert it by standard techniques 
order by order and end up with the series expansion of uo{ur) 

Uo{ur) =• + 



tins-* +2*-srr^ + i • • 



+ *-innr-» +3*-irTri~^ +2#-V • • A 
+ ^~rVT^ +^ • • • 4 +0«) 
We will further discuss the structure of this series in a later section. 



C. Resummed perturbation theory 

For two directed polymers interacting by a short range interaction (or equivalently one 
directed polymer interacting with a straight defect line by a short range interaction), an 
exact implicit equation for the dependence of the free energy per length from the interaction 
constant can be given by resummation of the perturbation series. 

In order to get these results, we review here the summation technique given in 02] and 



generalize it to arbitrary boundary conditions. After we got the implicit equation we will 
study some useful examples of specific boundary conditions. 
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1. Resummation for arbitrary boundary conditions 

The main idea, which leads to the summabihty, is that the coefficients of the perturbation 
series of the partition function have a product structure, which leads to a simple geometric 
series if they are properly decoupled. 

This decoupling is achieved by Laplace transforming the constituents of equation (p7D. 
We will call them for simplicity 

/(t) =Gi(0,0), M = JdriJ dr2GL|,(ri,r2), 



git) = I drGtir,0), and h{t) = J drGtiO,r). 

Their Laplace transforms are denoted by /, g and h respectively. 

Performing the Laplace transformation on the coefficients of the partition function, an 
implicit equation for the free energy per length can be extracted. The argumentation is 
given in appendix 0. The result is that 



lim EM = -zo. (33) 



-L| I — >oo 



where zq is the solution of 



l + gh^) = (34) 



with the largest (absolutely smallest) real part. 

This is an exact implicit equation for the free energy per length derived from perturbation 
theory. 

Introducing again dimension-less coupling constants we get the equation 

1 + UoL-J{L-\-fo - un)) = 0, (35) 

where /o denotes the dimension-less free energy per length of the free (mq = 0) problem. 
As it is obvious from the above equation, it can be calculated as the smallest pole of the 
function f{—L~^z). 

To further improve this equation, we write it as 

1 + uoF{ur) = (36) 

with F{z) = L~''/(L^^(— /o — z)). Since this equation has to be fulffiled for uq = and 
ur = 0, F has to behave like —1/z at zero (which can also be verified for specific infrared 
regularizations). This means that H{z) = F{z) + l/z is a regular function at zero. Expressed 
by H instead of F, the implicit equation for uq{ur) reads 

Ur-Uo + UqUrH{ur) = (37) 
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2. Results for specific infrared regularizations and the duality relation 



Now we want to calculate H for specific infrared regularizations, in order to compare the 
general equation (^) with results from the literature. 
Therefore we first reconsider the definition of H[z) 

roo fn+z 1 

H{z) = U / Gf(0, 0)e^*dt + - (38) 

JO z 
roo p 

= / Li-^G',i(0,0)e^«^- 



From dimensional analysis Gsl(0,0) must have the form 

G,L{0,0) = {sLy-'gin{sL). (39) 

Since Gsl(0,0) has to decay exponentially with the ground state energy of the quantum 
mechanical problem, i.e. like e"*-^", it is easy to see that H{z) is well-defined for 2; ^ in 
the infrared, because it follows that 

L'-'-'^GsiiO, 0)6^"' - 1 = s'-'hin{s) (40) 

with a for large arguments exponentially decaying function hjji{z). 

Also it is clear that ultraviolet divergences arise for e ^ 0. The derivatives of H with 
respect to z are less divergent in the ultraviolet regime; the k-th derivative will develop 
ultraviolet divergences at e ^ —k. Thus in the case e' ^ which we are interested in, also 
the first derivative of if at 2; = is ultraviolet divergent. 

The easiest case is the harmonic regularization, because the full propagator of the quan- 
tum mechanical harmonic oscillator is analytically known in all dimensions [^]. From this 



we extract the return probability for our directed polymer problem by inserting r = r' = 
and get 

L^~'Gsl{0, 0) = (2 sinh(7rs))~(^-^) (41) 



From the decay at large s we conclude /o = vr 



The Laplace transform of this function can be explicitly performed ( ||2^ 3.541) and gives 



27r r(e- 2^) 

For d = 1 or e = ^ this formula together with the equation (|36D coincides with the one- 



dimensional transfer matrix result in |2^. We want to stress that we just calculated the full 
transition function for the free energy from the Gaussian to the non-Gaussian fixed point 
in any dimension in the case of a harmonic boundary condition. 



It is remarkable that for = 3 or e = — 1/2 the equation ( |36D is equivalent to the d = 1 
case, if one replaces (7 up to a numerical factor by —l/g. This shows a duality (which exists 
for all dimensions which are symmetric with respect to c? = 2): The transition from the 
Gaussian fixed point in a dimension below d = 2 via an repulsive interaction to the non- 
Gaussian fixed point is exactly the same as the transition from the Gaussian fixed point in 
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the symmetric dimension above d = 2 to the non-Gaussian fixed point via an attractive 
potential. 

An exphcit calculation of the Laplace transform of the propagator with hard wall (Dirich- 
let) boundary conditions (see appendix |^) in (i = 3 gives the equation 



- = ^J2{un + fo)cot{^2{ut^ + fo)) (43) 



which shows the same duality relation with the transfer matrix result for d = 1 in |25 



D. Comparison of the two perturbative approaches 

We will now see that the resummed and the naive perturbation series after all manipu- 
lations presented in the last two sections are exactly equivalent. 

The resummation method gives us the exact relation (^) between the finite size coef- 
ficient ur of the free energy per length and the dimension- less coupling constant uq. If we 
assume that H can be expanded in a Taylor series around z = 0, a perturbation series can 
be easily extracted from this equation and it is the same perturbation series as we already 
calculated directly, if we identify 



1 



Jn— 1 



n intervals 



(n-l)\ dz 



n—l 



H{z). 





Using the specific representation ( |55D of H{z) by the function hi^^z) (pOj), which is directly 
connected to the return probability of the free problem, the coefficients are calculated as 

s^-^+^-^hiR{s)ds (44) 



in -I 



For the comparison with the transfer matrix results, we have to calculate the free energy 
per length at the transition point. Because we did use a renormalization point condition, 
which gives the renormalized coupling constant a physical meaning, in this framework, the 
finite size coefficient of the free energy per length is just the fixed point value of the renor- 
malized coupling constant. 

This fixed point value is the root of the (3 function, which describes the fiux of the 
renormalized coupling constant, and which according to the chain rule is given by 

/3M-^^. (45) 



Using equation (|36|) this can be expressed as 



This enables us to calculate this /5-function explicitly in the case of a harmonic infrared 
regularization, where the Laplace transform of the return probability can be performed (^2[). 
The result is 
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with the digamma function '^{z) = ^lnr(z). The behavior of the /3 function already has 
been shown in Fig. |^. For e < there is always one negative root and the slope of the f3 
function at the root is — e. It is remarkable that the (3 function has a well defined limit 
function for e —>■ —1 (e' ^ 0), which is the case we are mostly interested in. The limiting (3 
function is just the straight line uji + 27i. 




FIG. 2. The P function of a directed polymer in a harmonic potential. 



E. The limit e' ^ 

In the last section, we remarked that despite of the singularities in the coefficients of the 
perturbation series, a limiting f3 function for e' — exists. We heavily used the fact that 
we could calculate a closed form expression for the P function, and its properties are a bit 
awkward. Especially /5(0) is not zero any more in contrary to all finite orders of perturbation 
theory. 

If we want to get results for theories, where no exact solution is possible, we have to find 
a way, to extract the behavior of the jS function in the limit e' — *• directly from the lowest 
order terms of the perturbation expansion. 

We will give here first a general description of the regularization scheme and then show, 
how it can be applied to reproduce the limiting /3 function just shown. Moreover we will get 
results for other boundary conditions, too. 
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1. General regularization scheme 



Let us assume a perturbation series 



u ' 



(48) 



n=l 



the coefficients of which diverge in the hmit e' ^ 0. Assuming that e' = is a pole of hn of 
maximal ra-th order we use the Laurent expansion of 6„ 



(49) 



fe=0 



where lim^/^o Cra(e') = 0. We insert this expansion into the perturbation series and change 
the order of summation to 



fc=0 



.n=0 



oo 

+ E Cn(e)M^ 

n=l 



(50) 



If now the analytic continuations of the functions 



hk{x) = E 



(51) 



n=0 



have well defined limits for x ^ oo, the limiting form of f^i will be 



/o(^) = E 

fc=0 



lim /itfxl 



u 



(52) 



Moreover using just a finite number of these terms will give us a systematic expansion of f^i 
for small e'. 

To summarize the method, we have to extract just the most divergent parts of all the 
coefficients and use them to set up the functions hk- Those then have to be analytically 
continued and the limit a; ^ oo to be taken. The limits will give the coefficients of the 
limiting functions. 

In the following discussion this general approach has to be modified in some special cases, 
but the general idea will remain the same. 



2. The generic case of two directed polymers 

In both perturbative approaches, we have seen that the coefficients of the perturbation 
series can all be expressed by some set of constants oi, 02, . . . which can either be regarded 
as the fundamental graphical element of one arc in a diagram or as a derivative of the 
function H{z) defined above. The structure of the perturbation expansion is independent of 
the infrared regularization chosen. Different infrared regularizations just result in different 
values of the coefficients Oj. 

In the limit e' ^ only the ffist two coefficients ai and 02 will diverge, whereas all 



the higher coefficients remain regular, as it has been stated in section II C 2. At this point 
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we should remember that in the usual case e — only the coefficient oi is divergent. So 
the main difference in this problem is the existence of two divergent fundamental diagrams 
instead of one. 

To extract the pole structure, we assume the Laurent expansion of all the fundamental 
coefficients: 

dk = + akfi + ak,ie' + ak,2(^''^ + • . . , (53) 

where ak-i = for k > 3. 

Since we can calculate the perturbation series up to quite high orders, we can insert these 
expansions in the series and study the behavior of the most divergent terms. It shows up 
that the correct series to regularize is the series for , because this series has a pole of 
n-th order in e' among its n-th order diagrams. 

The singular parts happen to have a quite simple structure. All diagrams of a given order 
of divergence (order in ur minus order of the pole in e') can be expressed as a geometric series 
modified by coefficients polynomial in the running index n. The degree of those polynomials 
is at most the divergence order itself, as can be checked to quite high orders. So these series 
can be resummed, analytically continued and their limit for a; ^ oo can be taken. 

The explicit form of this series can be found in appendix 0. The regularized series in 
the limit e' = turns out to be a geometric series, that is resummed to 

hm ^ = ^ (54) 

e'^O e' ai _i + 02-1 Ur 

which leads to a linear limiting f3 function of 

lim (3{ur) = ^^i^ + Ur. (55) 

e'^O 02,-1 

The result is quite remarkable, because it states that the (3 function is linear independent 
of the infrared regularization chosen in the limit e' —>■ 0, which is surely not true for finite e'. 

The result of course reproduces exactly the explicit form of the /3-function, that we 
calculated for harmonic boundary conditions. 

It is interesting to state that the coefficients ai,_i and 02,-1 in the limiting function 
can be calculated much more easily than all the other coefficients, which were temporally 
involved in the calculation. We just have to recall the definition of ai and 02 (|4^). Integrating 
by parts the integral for ai two times and the integral for 02 one time and dropping the 
boundary terms, which are zero at least for e > 0, we get 

1 



Oi 



IWo 



00 



''h'iRis)ds (56) 



1 r°° I 

a, = -- KR{s)ds. (57) 

Since the divergences in e' are explicit in these formulas, we get the coefficients ai,_i and 
02,-1 by substituting e' = in the converging integrals. This results in 

ai,_i = h'j^{0) and 02,-1 = ^/_r(0). (58) 

So these coefficients are directly expressible by the return probability without any integra- 
tions. 
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3. The case of a vanishing leading singularity 



The above derivation has one severe problem. Although the calculated regularized series 
(and therefore the (3 function) is well defined in the limit of — >■ 0, we relied during the 
whole calculation heavily on the fact that ai _i is not zero. 

Unfortunately this is not true in the very important case of von-Neumann or periodic 
boundary conditions. These boundary conditions are the most simple ones and the pertur- 
bation series of more than two directed polymers can only be set up for periodic boundary 
conditions. 

With periodic boundary conditions the propagator is just the sum of free propagators 
connecting equivalent points 



1 \ 2 (r-r'-fcLj^)^ 

e~ 



r=0,r'=0 



Using L = \p2TiL\ as it is defined by (<l>(0))o = % we get (since /o = 0) 



At d = 4 /i/h(0) = 1, but /i;^(0) = 0. 

Since we know the whole perturbation series, we can just insert ai^_i = and repeat our 
above calculations. The structure of the poles changes, because now the only divergence 
resides in the arc, which spans 2 intervals. This leads to the fact that the order of the poles 
in e' increases only every second order and the coefficient connected to all poles is 02,-1- We 

2 

take care of this fact by grouping together all singularities with the same powers of 
The structure of the singular part of this series is then very similar to the previous one and 
the regularization scheme can be used in the same way as before. The explicit form of the 
important parts of the perturbation series is shown in appendix 0; in the limit e' — we 
get the expected result 

lim^^ = ^ + 0(4), (60) 

which is the same as if we had taken the limit ai _i — > in our above result. Therefore also 
the limiting /5 function is 

/?(«/?) = + 0(4) (61) 

There is one more remarkable thing in this (5 function: If we just add up the geometric 
series of the most singular terms and do not perform the limit e' 0, but just insert small 
e' we get the /5 functions shown in Fig. ^. From this we can extract two points. First, we 
can see, how the /3 function with a negative slope at zero approaches its limiting shape with 
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a slope of 1 at zero (and everywhere else). Moreover we remark that all three truncated 
f3 functions (they differ in the divergence order where they are truncated) have a common 
root at ±02 -i\/e'- The fact that this is true to all calculated orders strongly suggests that 
these roots are exact, which means that the fixed point approaches zero as 

UR ~ (62) 

which is just equation (p!0|). 



P(Uo) 




FIG. 3. Approximations to the /? function of a directed polymer with periodic boundary con- 
ditions for a small e'. The straight line is the limiting /? function (3{ur) = ur for e' = 0. The 
other three are calculated with the same e' = 0.005 but truncated at different orders. Two of them 
already coincide. 



4- Comparison between the perturbative and the transfer matrix approach 

For harmonic boundary conditions, the zero of the f3 function calculated by the regular- 
ized perturbative approaches is exactly equal to the finite size coefficient calculated in the 
transfer matrix picture. 

Since the von-Neumann and Dirichlet boundary conditions in the transfer matrix ap- 
proach and in the perturbation series are not absolutely equivalent (in the transfer matrix 
approach the boundary has a spherical shape), those results are not quantitatively compa- 
rable. Nevertheless qualitatively also those results are the same as far as the limit e' — > is 
concerned: 

In the case of hard wall (Dirichlet) boundary conditions, the perturbation series predicts 
a j3 function with a zero at a finite negative value. This is also the result of the transfer 
matrix calculation due to the finite value of the ground state energy of the free system. 
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With periodic (von-Neumann) boundary conditions, the perturbation theory and the 
transfer matrix approach both predict that C* approaches zero as e' ^ and they both 
result in the square root dependence (|1^) of e'. 

We can therefore conclude that the chosen method of regular izat ion of the perturbation 
series reproduces the exact results of the transfer matrix approach and therefore is a valid 
regularization scheme. 

III. ARBITRARY NUMBER OF DIRECTED POLYMERS 

Now we will develop a diagrammatic expansion of the theory of an arbitrary number of 
directed polymers N interacting by short range interactions. We will see that the structure 
shows strong analogies to the case of only two directed polymers, which suggests that the 
perturbation series can be regularized in the same way as the perturbation series for N = 2. 

A. Diagrammatic expansion 

1. Perturbation series of the partition function 

In order to keep the calculations as simple as possible, we choose periodic boundary 
conditions as the infrared regularization. That means that we identify the points r and 
r + kL± for every k G Z"^. Writing down the corresponding propagators, it is clear that 
the problem of N directed polymers with periodic boundary conditions is equivalent to the 
problem of N free directed polymers with periodic interactions. 

So the perturbation series for the partition function remains the same as p6| ) with the 
definitions 

*(^) = ( E E ^i^^it) - ^jit) - k^,LM (63) 

\i<j fci.eZ'* / 
and the free expectation values 

(r" r" Lii) r N 

r r r r ^'^^'-''^n'^W) _! r/-!! y ^2(t)dt 

(■)o = Ll^-^y dM ...J d% ydV/ . . .yd% J Vr,...Vr^.e - ' (64) 

As in the above discussion, we first have to calculate the multi-point functions in the per- 
turbation series of the partition function and combine them afterwards to the connected 
multi-point functions which are the integrands of the perturbation series of the free energy. 

Since the argument of a multi-point function is the sum over different 6 interactions, 
all sums can be extracted from the expectation value. In the n-th order the sums over 
i < j correspond to the different ways, the n interactions can be arranged among the 

directed polymers. Obviously a lot of arrangements are equivalent since for example 
interchanging the polymers that interact should not influence the expectation value of the 
sequence of interactions. So every possible arrangement of interactions will be accompanied 
by a combinatorial prefactor, which counts the number of equivalent arrangements. The 
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sum over all i < j can then be written as a sum over all possible arrangements with each 
arrangement multiplied by a convenient combinatorial prefactor. 

Since there are obviously possibilities to place the first interaction, all these 

prefactors will be multiples of This is the reason, why we divide the free energy by 

this factor in order to get the renormalized coupling constant. 

To improve the bookkeeping, we will represent every arrangement of interactions graphi- 
cally by drawing parallel lines representing the directed polymers and a connection between 
two of them for every 6 function between two polymers. The first connection belongs to 
the "time" ti, the last one to the time tn- To write down some simplified expressions later 
on, we also introduce a dotted line, which can represent a "time" in the diagram, where no 
interaction is present, and which is used to keep the relation between the time variables and 
the interactions, if an interaction is dropped due to a simplification rule. 

Between two "interaction times" we have to insert the free propagator for directed 
polymers which of course factorizes in a product of free 1-polymer-propagators 

»(r",At) 



GAtir',r")= f ' Pre-^/c/'^'Wdt 

J{r',0) 
1 



(27rAt) 



Since this is a Gaussian, all integrations over the can be performed. Moreover it is 
clear from the translational invariance of the propagator that only a dependence on "time" 
differences will arise. Therefore we change our integration variables. Instead of integrating 
over all ordered "times" < ti < ... < tn < we integrate over all existing "time" 
differences ti — In the limit of infinitely long directed polymers (L|| —>■ oo) the domain 
of integration is from zero to infinity for each of these variables, whereas the integration over 
ti just gives a factor of L\\ which is canceled because we want to calculate the free energy 
per length and per pair of directed polymers. The new integration variables will be made 
dimension-less by a factor of L and we call them si, . . . , 

To calculate the integrands, some simplification techniques can be used and especially 
two general simplification rules, which we will call lemma A and B can be derived. The 
integrands of the first three orders of the perturbation series for the partition function can 
be explicitly calculated by using these techniques. The detailed simplification rules and the 
calculation for the first three orders can be found in appendix |y and lead to 

-{m)o = L~/, (65) 



A^(A^- 1) 

and 



N{N - 1) 

($(0)<l>(Lsi))o = L-/RPis,) + ^ ^^ ^LI^^ (66) 



^-Ll"'R^{s2) (67) 



(iV + l)(iV-2)^_2d: 
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{N - 2){N - 3)(iV^ + 3A^ + 4) ^ 
4 



-3d 



with the two-polymer return probabihty 




(68) 



We want to stress at this point that up to the third order the only quantity arising is the two- 
polymer return probability -R^(s) and the only effect of more than two directed polymers are 
the combinatorial prefactors. Unfortunately we cannot just stop our discussion here, because 
in the forth order, generically new diagrams arise. Since the number of diagrams heavily 
increases, it is convenient to use some computer algebra in order to generate all diagrams. 
In the sixth order for example there are 29388 possible arrangements of the interactions, 
which can be reduced to 5300 by automatically applying lemma A. 
The new diagrams in the forth order are the following ones 



I 



T 



They can of course also be integrated out, which gives L^'^KF{si + 2s2 + S3) for the first 
one. The other two contain two sums over Z*^ that cannot be decoupled. The second one 
for example reads 



det + 



(kk') 



Jti + t2 + h ¥2 



t2 



4det| 



oS2 



(69) 



Here it is already written in a form that makes it obvious that it is connected to a two- 
dimensional theta function and that it can be simplified by applying the corresponding 



transformation formula to 



1 

Z2 




Calculating the diagrams to higher orders shows us that there is a simple recipe, that allows 
to perform all spatial integrations formally, if the form of the diagram is given. The rules 
are as follows: 
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1. 



Mark all loops in the diagram that are necessary to pass each interaction at least with 
one loop. Assign an orientation to every loop. 



2. If there are m loops necessary, use an m x m matrix T and identify each row and 
column with one of the loops 

3. pass for all "time" intervals i over all lines of the diagram and 

• add |si to every diagonal element of a loop, that passes the line 

• add ±|<Si to both of the off diagonal elements of two loops, that pass together 
through the line. If they pass in the same orientation the plus sign has to be 
used, otherwise the minus sign. 

4. The value of the diagram is then 



order 

L 2 



(70) 



This prescription also implies the vahdity of the lemmas A and B used above. 

Since there are several possibilities to choose the loops in a diagram, the correspondence 
from a diagram to a T-matrix is not unique. For example in the forth order diagram, the 
loops could be chosen as: 







T 


( 


— 












3 





which leads to the T-matrices 

( h+h \t2 \ . ( h+t2+h lt2 + h\ 

[ 1^2 h + ts) [ lt2 + ts h+h )■ 

Although the T-matrices are not identical, the sums, which determine the value of the T- 
matrix are the same, because the second of those matrices is reproduced, if in the sum with 
the first matrix, the summation variable k2 is shifted by ki. During all those equivalence 
transformations, the determinant of the T-matrix docs not change, so that two T-matrices 
from different diagrams can only be equivalent, if their determinants are identical. Unfortu- 
nately there exist also T-matrices which arc not equivalent but have the same determinant. 

Since the whole procedure described above is purely mechanical, it can be implemented 
by computer algebra. The only thing left to do is manually checking the equivalence of T- 
matrices which have the same determinant. This reduces the number of different diagrams 
in the fifth order from 348 to 88. 
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2. Perturbation series of the free energy 



Since we now know how to compute the multi-point functions, we can calculate the free 
energy using the formula from appendix p. Analogous to the connected two point functions 
we used for the 2 polymer problem, we have to express everything in terms of correlation 
functions, that decay properly for increasing time differences. In analogy to the definition 
of Qi, which in our terminology reads gi{s) = L'^Rp{s) — 1, we will use functions of T- 
matrices, where all sub-matrices of lower dimensions are subtracted with alternating signs, 
as for example 



g2{a,b,c) = 












f-l 




[fcez 




fcez J 



and analogously for a higher number of loops. 

Of course the subtracted terms depend on the T- matrix itself and are different, if one 
chooses different T-matrices which represent the same diagram (which give the same value for 
the first term of the above sum). If one chooses the wrong representation, one will discover 
that there are terms, the time integrals of which will not converge. But it is always possible 
to choose a representation, that leads to terms which are each integrable (of course the sum 
of all terms is integrable in every case, it is just possible that one chooses an inconvenient 
partition of the integrand.) 

Once the coefficients of the free energy are represented as a sum over time integrals, 
further simplifications can be performed, by multiplying the integration variables by constant 
factors and applying all kinds of relations as the ones discussed in appendix y. 

We will continue to represent every of those integrals by one of the diagrams, which is 
responsible for the main term of the integrand. Since the integrands themselves depend on 
the way, the loops are chosen, we must specify our choice. Those parts of diagrams, which 
can be expressed only by gi (i.e. which are pure 2-polymer-diagrams) will be represented 
by the same diagrams as they have been used in the 2 polymer case. 

The resulting diagrammatic expansion for the free energy is shown in Fig. ^ up to the 
fifth order. 
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FIG. 4. Diagrammatic series expansion of the free energy per length up to the fifth order for 
an arbitrary number N of interacting polymers. 



It should be stressed that the number of diagrams involved in this final expansion is 
very small compared to the number of diagrams in the original perturbation series of the 
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partition function. Since the applicability of most of the simplifying formulas as those from 
appendix relies on very specific relations between the coefficients of the diagrams involved, 
it is very unlikely that this structure evolves just by chance. This is a strong hint that there 
exists a similar equation as (^) for an arbitrary number of directed polymers, that produces 
this relatively simple structure of the perturbation series. Unfortunately no such formula 
could be found up to now. 



B. Comparison with the N = 2 case 

We can now compare the perturbation series to the perturbation series of two directed 
polymers. 

First there are a lot of diagrams which appear already in the series of two directed 
polymers, that now have prefactors with a polynomial dependence. We will call them 
"linear" diagrams, because they have no nested loops. The value of one arc itself is of course 
the same as for two directed polymers. Especially only the arc, that spans two intervals has 
a pole at e' — > 0. This pole is of the first order. 

Since we have explicit expressions for all of the diagrams, it is also possible to extract 
at least the order of the pole of each of the nested diagrams as e' ^ 0. This can be done 
by extracting the UV dependence of the integrands and partial integration. It comes out 
that all nested diagrams are less divergent than the product of two-interval arcs in the same 
order of the perturbation theory. 

Therefore, as in the case of two directed polymers, the powers of the two-interval arc are 
the most divergent diagrams and the order of the pole grows by one only in every second 
order of the perturbation series. 

Although we gave here the diagrammatic expansion of the series ur{uo), the same state- 
ments are true for the inverted series uq{ur), since it consist of exactly the same diagrams 
with other prefactors. 

Analogous to the case of two directed polymers, we write the partial sum over all those 
most divergent diagrams as 

°^ = — fren I — ^ I + l^ss divergent terms (71) 

where freg is an even function, with Taylor coefficients, that stay regular as e' — >• 0. If we 
now assume that the fact that the series for freg{x) can be analytically continued to x — cxo 
and that lima;^oo freg{x) exists is not a particularity of 2 directed polymers, we get a leading 
behavior of 

^ ^ with a = hm U,{x) (72) 

e Ur x— >oo 

Differentiation then leads to the limiting (3 function 

P{ur) = Ur + higher order terms. (73) 

The slope of 1 of this (3 function shows that = is not just the Gaussian fixed point, 
which would have a slope of e = — 1. Thus we expect that for e' ^ the non-Gaussian fixed 
point goes to zero for each N. Since the /? function is a regular function of UR/^/e' also the 
behavior proportional to is independent of A^. 
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C. Consequences for the limit ^ 

The independence of the behavior of the fixed point near d = 4 renders the hmit N ^ 0, 
which describes the directed polymer in a random medium, trivial Thus we easily conclude 
that the finite size amplitude of the free energy per unit "time" also vanishes as given in 
equation ([T0|) for the directed polymer in a random medium. 

This shows that the singularities in the perturbation series, which arise at d = 4 are 
not only a formal problem of the approach, but they lead to a non-analytic behavior of a 
physical quantity. 

Although we are surely not able to access the strong coupling fixed point of the KPZ 
equation with these perturbative methods, the shown vanishing of the weak coupling fixed 
point, which is under the Hopf-Cole transformation equivalent to the studied interaction 
fixed point in the directed polymer picture, is a strong hint that 4 is really the upper critical 
dimension for the KPZ problem. 

Although the duality relations between dimensions d and 4 — ci in section [11 C 2| only have 
been observed in the case N = 2, one could speculate that such a relationship still exists for 
arbitrary and the KPZ problem. This is especially interesting since Frey and Tauber find 
in p9| , p!0[| that their coupling constant approaches zero in the limit d ^ 0. Surely it is not 
absolutely clear how the two coupling constants are related and this point deserves further 
investigations. 

IV. CONCLUSION AND OUTLOOK 

We have studied the problem of directed polymers with short range interactions. The 
main interest was the behavior of the directed polymer system in the limit — > 4, where the 
normal renormalization treatment, which expands around d = 2 breaks down. Since for two 
directed polymers exact transfer matrix results are available, we developed a regularization 
procedure for the perturbation series in the limit d A in the special case of two directed 
polymers. We could prove that there is an agreement between the transfer matrix results and 
the predictions of the regularized perturbation series for two directed polymers. Moreover 
an exact implicit equation which generates the whole perturbation series and its connections 
to the diagrammatic expansion has been studied. 

For an arbitrary number of directed polymers, a diagrammatic expansion of the free 
energy up to the fifth order has been established. Although it stems from a large number 
of terms this expansion is relatively simple. This suggests that it can also be generated 
by a simple equation, which has, however, eluded us so far. Since the pole structure of 
the perturbation series is very similar to the pole structure in the case N = 2, the leading 
behavior of the (3 function in the limit d —>■ 4 has been derived. It shows that the interaction 
fixed point approaches zero as d ^ A proportional to independent of A^. Therefore 
the weak coupling fixed point of the KPZ equation (which is mapped to the problem of a 
directed polymer in a random medium) is also expected to approach proportional to y/e' 
in the limit c? — 4. 

For d > 4 the transfer matrix calculations show for two directed polymers, that the 
finite size amplitude of the free energy per unit "time" at the unbinding transition develops 
anomalous scaling behavior analogous to that of the free energy in the bound state of the 
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infinite system |1T9| , |20[] . The finite size amplitude therefore depends exphcitly on the short- 
distance cutoff. Thus it is no longer possible to define a universal quantity from the free 
energy. If this behavior persists for arbitrary A^, the critical behavior of the KPZ equation 
at the roughening transition above d = 4 is less universal than below 4 dimensions. 

We conclude by noting that the unbinding transition of semi-flexible polymers (i.e. di- 
rected lines which are governed by the curvature instead of a line tension) with a short-ranged 
attraction of parallel polymer segments in one transversal dimension formally can be mapped 
on the problem of two strings in four transversal dimensions. Therefore the method studied 
here can also be applied to examine this unbinding transition 
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APPENDIX A: SOLUTIONS OF THE RADIAL SCHROEDINGER EQUATION 

WITH CONSTANT POTENTIAL 

We have to calculate the asymptotic behavior of the ground state energy, which belongs 
to the rescaled version 

r (d-3)(d-l) 1 ^ , 

+ ^ — — + viy)^m = Eom 

of the radial equation ([M])- The rescaled potential is 

where 0(r) = $(ra), Eq = a^Eo, V{r) = a^V^(ra) and Vq = o^Vq. 

Obviously the wave function must consist of general solutions of the Schroedinger equa- 
tion with constant potential. 



The general solution of this differential equation is 



<f){y) = I Ay'^ + By^ E = ■ 

[ A^I^^yV^) + B^K\,\{y^) E < 

Because the wave function should be regular at the origin at least for d > 2 only the first 
solution is possible for the interaction region. This means that for y < 1 

( A^J\,\{yVEo + Vo) Eo + Vo>0 
<P{y) = A^y\^\ + K) = . (A3) 



1 



^v^/h(W-^o-^o) ^0 + K) < 
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In the outer region the boundary conditions at y — L/a lead to the wave functions 



0(y) = B[^J\4y^Eo) - gi4^Eo-)^N\4y^Eo)] 

a 



for Eq > 0, 



- (^^j"* %V] Dirichlet 

von- Neumann 



By^-i^ 



for £'0 = 0, and 



for £'0 < respectively with 



f\e\{x) 



Dirichlet 
von- Neumann 



and 



9\e\{x) = < 



_J\e\ + lix)_ 

C Dirichlet 

if'^'"'"^^^! von- Neumann 



(A4) 



(A5) 



The total solution consists of these two solutions if they obey the condition of continuous 
differentiability, which is most conveniently written as 



0(1-) = 0(1+) 



and 



1, 



</>'(!-)- (kl + 2^1-) 



^\l+)-i\e\ + -mi+). 



Since this are two conditions for the two free coefficients A and B, we only get a solution 
for energies, where the determinant of the linear system of equation vanishes. This gives us 
the equation connecting Eq with Vq and 
This equation will always have the form 



±En 



(A6) 



with some functions r and s. Since the ground state energy is expected to vanish for 
L± — 00 at least proportional to {a/L±y, we have to distinguish two cases which are 
determined by the value of Sq = ^^t^x^os{x). If So is a finite value, Eq will asymptotically 
decay as ±c^{a/L±y, where c is the smallest positive solution of r{x) = Sq. So is exactly 
the coefficient, we are interested in. If sq = 0, the constant c is a root of r{x). Since r(0) = 0, 
c can be zero, but only if the signs of r and s for small positive arguments are equal. If 
not, c is the smallest positive root of r. If c is zero, one can expand r for small arguments 
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and furthermore extract the leading behavior in a/L^ of £'0, which then decays faster than 
quadratic. 

Since we are especially interested in the asymptotic behavior of the ground state energy 
at the phase transition point K, we have to identify it. It is defined by the condition 
that the ground state energy in an infinite system approaches zero from below. Since we 
know the solutions of the Schroedinger equation in the interaction region and outside (for 
Lj_ — >■ cx) only the regular solution proportional to K\^\ is possible), we can exploit the 
matching conditions for + K > and Eq <Q and take the limit £^0 ~^ 0, which gives 



We now can systematically apply the above scheme to all possible combinations Vq = 0, 
< I/q < K and Vq = K, £0 < -Vq, Eq = Vq, -Vq < Eq <0, Eo = and Eq > and both 
boundary conditions and extract all possible values of the coefficient c. The smallest one is 
the ground state energy. 

The only point, which has to be handled with care during this calculation is the series 
expansion of the different Bessel functions involved. Since especially for Vq = K the leading 
terms of the expansions cancel, the sub-leading terms have to be used. But the sub-leading 
terms arc of totally different origin if |e| > 1 or |e| < 1. This produces the difference of the 
ground state energies in d > 4 and d < 4 as we expect them. 



Since we need it during the calculations, we give here the expansion of the free energy 
per length written by n point functions up to the forth order. It is: 




2|e|. 



(A7) 



APPENDIX B: EXPANSION OF THE FREE ENERGY 




Fjg) - m 

L|| 



gmo))o-g 



[■00 

' {mm)o - {m)idt + 

Jo 




il 



dtidta ($(0)$(ti)$(t2))o - (^'(0))o($(ti)$(t2))o 



0<tl<t2 



($(ti))o($(0)$(t2))0 - ($(t2))o(*(0)$(ti))o + 2($(0))3 







0<tl<t2<t3 



-($(0)$(i2))o($(ti)^(i3))o - ($(0)$(i3))oWi)^(i2))o 

-($(0))o($(tl)$(t2)$(i3))o - ($(0))o($(0)$(i2)$(t3))o 
-($(0))o($(0)$(ii)$(i3)>0 - ($(0))o($(0)$(ii)$(i2))o 

+2($(0)$(ii))o($(0))g + 2($(0)$(i2))o WO))g 



-r2{^{t,)<^{h)Um)l + 2{^{h)^{h)Um)l - 6($(0))^ 
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In the case of only two directed polymers, this whole series can be reexpressed only 
by connected two point functions and then represented diagrammatically as explained in 
section 111 B. 



It reads then: 
mr(mo) 



APPENDIX C: PROVING RELATIONS AMONG DIAGRAMS 

Every diagram is related to a multiple "time" integral over a product of connected two 
point functions, the arguments of which are sums of the different integration variables. If 
we assume that the two point function can be Laplace transformed, general rules among the 
diagrams can be proved. We will show with one example how this proceeds. We introduce 
the Laplace transform gi of gi. Then: 

/■OO roo poo roo roc poo 

/ dsi ds2 ds3gi{si + S2 + S3)gi{s2)+ dsi ds2 ds3gi{si + S2)gi{s2 + S3) = 
Jo Jo Jo Jo Jo Jo 



{2Tci) 



dzi / dz2gi{zi)gi{z2) / dsi 



00 roo 

dS' 



ds-^ 



2 2 



+ e 



^l(si+S2)+22(s2+S3) 



From its definition c is a path parallel to the imaginary axis with positive real part. In 
order to be able to interchange the integrations, we have used that the inner integrals exist, 
which is of course only true, if the real parts of the Zi are negative. But since the connected 
correlation functions decay exponentially for large arguments, the Laplace transform gi is 
still analytic in some region to the left of the imaginary axis and the integration contour can 
be shifted there without changing the value of the integral. 

After doing so, we can evaluate the integrals over the exponentials and get 



(27rz)2 ' ' dz2gi{zi)giiz2) 



dzi / dz2gi{zi)gi{z2) 



(27r2)2 
1 



2zl{Zi+Z2) 2zl{Zi + Z2) ZiZ2{Zi + Z2 



dzi / dz2gi{zi)gi{z2) / dsi 1 us2 

Jo Jo JO 



c Jc 

00 /"OO 



(27r2)2 

roo roo roo 

/ dsi ds2 ds3gi{si)gi{s2 + S3 
Jo Jo Jo 



dso 



^g^^Z^(s2+Ss)+Z2 



^2, 



The diagrammatic expression for this equation is the one in section |IIB . 
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APPENDIX D: DERIVATION OF THE EXACT IMPLICIT EQUATION FOR 

THE FREE ENERGY 



Using the abbreviations from section [11 C| the n-th order term of the partition function 
series is expressed by 



/ dtn / dt„„i . . . / dti 
Jo Jo Jo 

9{h)f{h - h)--- f{tn - tn-i)h{L\\ - t„) 



Laplace transforming this with respect to L\\ yields g{z)f{z)"'^^h{z) for the Laplace trans- 
forms. 

Obviously the Laplace transformed perturbation series is just a geometric series and can 
therefore be resummed. After back-transformation we end up with 

A _ 1 = 9_ f IMM^.L,,.^^ 

Zo 2mU Jcl + gf\z) ' ^ ^ 

where c is a path in the complex plane parallel to the imaginary axis. Since we can obviously 
close this path by a circle at z ^ — oo, the integral is given as the sum of the residues of the 
integrand in the half plane of negative real parts. From the form of the integrand it is clear 
that all residues will be some prefactor times an exponential with the position of the pole 
times L|| as its argument. In the limit of L\\ — >■ oo only the pole with the smallest decay 
rate (i.e. the one with the smallest absolute value of its real part) survives. 

By construction, it is clear that TV has a leading dependence on L\\ of e~^ii^'^'^o, whereas 
the poles of g and h are exactly the negative eigenvalues of the Schroedinger operator corre- 
sponding to the free directed polymer problem (described by T^o)- The smallest eigenvalue is 
the leading term of In Zq itself, the contribution of which to the integral must cancel against 
the 1 on the left hand side. 

From that we conclude that the leading term of Z for large L\\ is some prefactor times 
exp(zo-^ll), where zq is the solution of the equation 

l+gf{z)=Q (D2) 

with the largest (absolutely smallest) real part. This decay rate is therefore the leading 
contribution to the free energy per length in the limit Lm oo. 



APPENDIX E: HARD WALL RETURN PROBABILITY 

To calculate the return probability of a l+ci- dimensional directed polymer in a round box, 
we can use the "quantum mechanical" expression of the propagator by the eigenfunctions 
of the "particle in a box" problem. The eigenfunctions of the particle in a box are Bessel 
functions of the first kind and for > 2 we get with the correct normalization conditions 



Gi(r,r') 



oo 



^± 1=0 Jrf_2 



E 



2-d 

r 2 Jd_2 



+1 



ar]r 



2 Ji^^i[ar' 



+1 



J d- 



'aLj 



(El) 
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which for c? = 3 is the heat equation kernel in p8 . 

In the hmit r 0, which we need for the return probabihty, only the / = terms 
stay finite. The sum over the for / = radially symmetric eigenfunctions of the angular 
momentum operator is just one over the surface of the d-dimensional unit sphere. Thus the 
return probability is 



r(f) 



2 



TT' 



-L J 



E 

<i-2 {ctL^ 



d-2 



)=0 



Jd\ 

2 



aL I 



(E2) 



Its Laplace transform can formally be calculated term by term, but since we know that 
it is ultraviolet divergent for ci > 2 we introduce a lower cutoff a of the integration, which 
give 



a 



d-2 



7r^/2(2L^)'^-2 



Jd-2(")=0 



Jd 

2 



a] 



(E3) 



In the prefactor the limit a — > is possible without any difficulties. 

If we now specialize to the case ci = 3, where the roots of the Bessel function are just the 
integer multiples of vr, we can insert the especially simple expressions for the Bessel functions 
and end up with 



oo 



2 2 
6 2 



n- 



=1 1 + 



2zL^, 




k=0 ^ " ' 71=1 

where we have absorbed a factor of 1/ L\ into a and omitted the geometrical prefactor. For 
> 1 the limit a — is possible and we get 

k 

'2 \k 



oo 



n=l 



oo 

-E 

n=l 



oo 

-E 

k=l 



2 2 

e 2 



TC 



The last equation can be verified by representing the values of the C function by Bernoulli 
numbers which gives exactly the series expansion of the cotangent function. The k = 
sum obviously diverges for a — >■ 0. If we add ^, it is half of the value of theta function 
at zero, which diverges like with no sub-leading algebraic terms. So, if we ignore the 
divergence, the first sum contributes — f to the return probability. Thus the regularized 
Laplace transform of the return probability is 



1 



(E4) 



if we add all the geometrical prefactors again. This leads to the equation for the free energy 
per length of the system with a short range interaction cited in the main text. 
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APPENDIX F: EXPLICIT REGULARIZATION OF THE PERTURBATION 

SERIES 



If we insert the Laurent expansions of the coefficients into the perturbation expansion 
of ""^J'"^ we get a regular part of this series to the leading orders in ur of 

ai,iu% + (02,1 + 2ai _iai,2 + 2ai,oai,i)M| + 0{u%) 

As discussed in the main text, the singular parts has a quite simple structure. All terms of 
the same divergence order can be combined to geometric series with polynomial prefactors. 
The most singular terms (the first three divergence orders) are 



+ 



ai-i 



n h (n — Ij- 



^{n + l)ai,; 



+ in 



+{n^ - n) 



^2,0 , Ct2,0 , 

n h n h 

ai -1 

— 3n - 



Ur + 
n{n + l) (ai,o)^ 



Ol,0O2,-l 



Ol,-l 



+ 



2 (a2, 



'1 



u 



R 



0(4)} 



To complete our program, we just have to find the analytic continuations of series of the 
form 



E 

n=l 



(Fl) 



and their limit for a; — > cxo. This is easy because they all are derivatives of geometric series. 
It turns out that the limit for a; —> 00 is —1 for A; = and for all A; > 0. So we get the 
contributions of the singular terms in the limit e' — > by just inserting n = in above 
expression and taking the negative value of it. If we do that, we arrive at 



lim 



Uo{ur) 



1 ^ 02,-1 

+ 1 ^Ur 



(F2) 



K-l)' 2 , (Q2,-i) ' 

(ai,-i)^''^+(ai,-i) 



4^R 



0(4) 



(for the third order coefficient we need one term more in the above formula for the singular 
parts, which has been omitted, because it is easily computed but quite lengthy.) 

This is obviously the beginning of a pure geometric series. 

If Qi-i — 0, the regular part of the series UQ{uR)/e' is just 



di luj} + 0(u 



R) 



(F3) 



The singular part consist again of geometric series with polynomial coefficients and explicitly 

reads 
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0-2, 



-lUR ^1 V J 



"2-1 „=i 



n{n + 1) 



(ai,o)^ + 72 02,0 



+ n[n + l)ai,oa2,o + (''^ + l)ai,ia2 -i + '^o-s.o H 7 («i,o) 



6 



+ 0(4) 



The limit e' ^ is again performed by inserting n = and taking the negative value which 
reproduces the expected result (|60D. 



APPENDIX G: FIRST THREE ORDERS OF THE PARTITION FUNCTION FOR 
AN ARBITRARY NUMBER OF DIRECTED POLYMERS 



During the calculation of the integrands in the series expansion of the partition function, 
most of the terms can be strongly simplified by 

• using the symmetry of the one particle propagator 

• moving parts of the arguments of the one particle propagator from one argument to 
the other using the fact that the propagator depends only on the difference of the 
arguments 

• translating R'^-integrations by fcL^ terms 

• translating TJ^ summations by k' from other sums 

• combining of sums over TJ^ and integrals over [0, L\\^ to integrals over R'^. 

With this technique, it can be generally shown that a directed polymer which is not 
involved in any of the interactions does not contribute to the value of a diagram and that a 
directed polymer that is involved only in one interaction contributes just factor of . We 
will call this lemma A and represent it graphically as 




Moreover it is possible to prove lemma B: 




With this preparation, it is easily possible to compute the 1-, 2- and 3-point-function. The 
1-point-function has just one diagram with the prefactor A^(A^ — l)/2. 
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which has the value L^^ according to lemma A. Therefore 

^^i- ,G1, 



In the second order there are three types of diagrams with combinatorial prefactors of 1, 
2(N-2) and (A^ - 2)(A^ - 3)/2 respectively (omitting the general prefactor of A^(A^ - l)/2). 



X 



The last two are reduced by lemma A to L^^'^, whereas the first one has the value 



feezd 



Integrating over r results in 



d 



e 



where the second equation comes from the fact that the sum is the value of a theta function 
at zero |27]. 

Combining everything, we get in the second order equation (|66|) . The third order consists 
of 16 different diagrams. All of them but one can be evaluated by applying lemmas A and 
B and in the end we get equation 
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